Intermittency exponents and energy spectrum of the Burgers and KPZ equations 

with correlated noise 



Mahendra K. Verma * 
Department of Physics, Indian Institute of Technology, 
Kanpur - 208016, INDIA 
(February 5, 2008) 

We numerically calculate the energy spectrum, intermittency exponents, and probability density 
P(u') of the one-dimensional Burgers and KPZ equations with correlated noise. We have used 
pseudo-spectral method for our analysis. When a of the noise variance of the Burgers equation 
(variance oc k~ 2a ) exceeds 3/2, large shocks appear in the velocity profile leading to < |w(fe)| 2 >oc 
k~ 2 , and structure function < \u(x + r,t) — u(x,t)\ q >tx r suggesting that the Burgers equation 
is intermittent for this range of a. For — 1 < a < 0, the profile is dominated by noise, and 
0^ ' the spectrum < |/i(fc)| 2 > of the corresponding KPZ equation is in close agreement with Medina 

et al.'s renormalization group predictions. In the intermediate range < a < 3/2, both noise 
and well-developed shocks are seen, consequently the exponents slowly vary from RG regime to 
a shock-dominated regime. The probability density P(h) and P(u) are gaussian for all a, while 
P(it') is gaussian for a — —1, but steadily becomes nongaussian for larger a; for negative u' , 
f"*^ ' P( u ') x exp(-ax) for a = 0, and approximately oc u'~ 5 ^ 2 for a > 1/2. We have also calculated the 

energy cascade rates for all a and found a constant flux for all a > 1/2. 
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I. INTRODUCTION 



Stochastic Burgers equation has been studied extensively because of its close connection with the Navier-Stokes 
equation. Recently it has been found that Burgers equation with correlated noise shows intermittency [0-0], however, 
these calculations have been done for specific degrees of noise correlation. In this paper we vary the applied noise from 
uncorrelated regime to strongly correlated regime and study the energy spectrum of h and it, as well as intermittency 
exponents and the probability densities of u and du/ dx, where u(x, t) is the velocity function appearing in the Burgers 
equation, and h(x,t) is the surface height of KPZ equation. 
The one-dimensional Burgers equation fl is 



k> \ du , du d 2 u 

m +Xu dx = ^ + c > (L1) 



where u is the velocity field, A is the strength of the nonlinear term (A = 1 in the standard equation), v is the viscosity, 
and £ is the noise. We assume that the noise C(k,t) is gaussian (C(k) is the Fourier transform of ((x)), and in Fourier 
space follows distribution 

< C(fc, t)((k', t') >= 2Dk- 2a (2tt) 2 S(k + k')S(t - t') (1.2) 

Tatsumi and Kida ||], Kida j6|, Gotoh [Q, and Bouchaud et al. Q have solved the noiseless Burgers equation (( = 0) 
exactly in limit of large time t and zero viscosity. The solution is a linear profile with sharp discontinuities, which 
are called shocks. This shock solution yields structure function < \u(x + r,t) — u(x,t)\ q >cx r indicating that the u 
is intermittent. The solution also yields the energy spectrum < \u(k)\ 2 /2 >oc k~ 2 (u(k) is the Fourier transform of 

The noisy Burgers equation as well as the noisy Kardar-Parisi-Zhang (KPZ) equation have been studied by many 
authors. The Kardar-Parisi-Zhang (KPZ) equation, which describes a generic set of surface growth phenomena, 
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is closely related to the Burgers equation. The replacement of u by —dh/dx in the Burgers equation yields KPZ 
equation: 



dh A fdhy d 2 h „ _ . . 

where h(x, t) is the height of the surface profile at position x and at time t, A is the strength of the nonlinearity, v 
is the diffusion coefficient, and / is the forcing function. We again assume that the noise f(k, t) is gaussian, and in 
Fourier space follows distribution 

< f(k, t)f(k\ t') >= 2Dk- 2p (2vr) 2 S(k + k')S(t - t'). (1 .4) 

It is easy to see that |w(fc)| 2 = k 2 \h(k)\ 2 and a = p— 1. In the following several paragraphs we will describe the work 
done by various researchers regarding energy spectrum, intermittency exponents, and probability densities of h and 
u when noi se is present in the system. 



Eq. (1.3) has been solved by Kardar-Parisi-Zhang |9|] for a = — 1 and by Medina et al. [|l(J for —1 < a < using 
renormalization group analysis. They have calculated the roughening exponents xkpz and (3kpz, which characterize 
the dynamical properties of the equation, e.g., < |/i(fc)| 2 >oc k~ 2xKPZ ~ 1 . The RG analysis is not applicable beyond 
(7 = 0. Using the connection of the Burgers equation with the KPZ equation, the dynamical exponents of the Burgers 
equation is also automatically solved for this set of a. Recently, Chattopadhyay and Bhattacharjee (jflj have applied 
mode coupling scheme and obtained the same exponents for KPZ equation for — 1 < a < 0. 

For a = —1, Barabasi and Stanley |l2| have shown that in ID KPZ equation, the function dh/dx{= —u) has a 
gaussian probability density under steady-state. Since h! is gaussian, we can immediately determine the structure 
function T q (r) =< \h(x+r) — h(x)\ q > using the following arguments. Since Ah w rh' for small r, Ah[= h(x+r) — h(x)] 
will also have a gaussian probability density, i.e., 

Pm * expf-^Y (1.5) 

Therefore, 



/OC 
P{Ah)\Ah\ q dAh 
-OO 



(1.6) 



In the region where < |A/i| 9 > is a powerlaw (region of our interest), ay will also be a powerlaw, say, ay cx r x . 
Therefore, < \Ah\ q >cx r xq . 

The above discussion clearly shows that if the derivative of the function has a gaussian probability density, then 
the exponent of the q— th order structure function (denoted by (^ q ) is proportional to q. These kind of functions are 
called non-intermittent functions. On the other hand, if the probability density of the derivative is nongaussian, or 
the exponent of q— th order structure function is not proportional to q, then the function is said to be intermittent. 
The intermittent functions have typically powerlaw tails instead of gaussian tails in their probability densities, hence, 
probability of occurrence of large value of the function is higher for the intermittent functions then the non-intermittent 
functions. In other words, large events, which are absent in non-intermittent systems, occur in intermittent systems 
in a bursty manner. In fluid turbulence, it has been found that P(u) is gaussian, whereas P{u' L ) and P(u' T ) (where L 
and T denote longitudinal and transverse components) have powerlaw tails for small r. The exponents of the structure 
functions are not proportional to q either. Hence, fluid turbulence is said to have intermittency. These observations 
are explained using localized vortices. For further discussions on intermittency in fluid turbulence, refer to Frisch |13[1 
and references therein. 

Chekhlov and Yakhot , and Hayot and Jayaprakash || numerically solved the structure function of the Burgers 
equation with correlated noise for a range of a. For a = 1/2, Chekhlov and Yakhot |Q obtained < |AM(r)| 9 >oc r 
and claimed that the Burgers equation is intermittent according to the above definitions. Hayot and Jayprakash || 
found similar results for < a < 1/2. This behaviour was attributed to the shocks present in the system |l]||. 

Regarding the energy spectrum, Chekhlov and Yakhot pj obtained Kolmogorov's energy spectrum, i.e., |u(fc)| 2 oc 
fc~ 5 / 3 for a = 1/2. They argued the constancy of energy flux to be the reason for Kolmogorov's spectrum in the 
noisy Burgers equation with a = 1/2. In this paper we will show that energy flux is constant for all a > 1/2, yet 
the spectral exponent varies from 5/3 to 2. Hence, Chekhlov and Yakhot JjJ claim that the constancy of flux implies 
Kolmogorov's spectrum for Burgers equation is incorrect. 
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Hayot and Jayaprakash || have numerically calculated the energy spectrum for —1 < a < 1/2. For —1 < a < 0, they 
observed that at large k the behaviour is that of free field, and at low k is in agreement with the result derived using 
RG treatment. In our simulation we do not find any such crossover. For < a < 1/2, Hayot and Jayaprakash argue 
that the shocks determine the exponents; their numerical exponent for energy is quite close to the renormalization 
group formula even though RG is not expected to work in this regime. They also find that the exponent f^Burg = 1/2 
for the Burgers equation. 

Polyakov Q has applied methods of quantum field theory to calculate probability density P{u') of the Burgers 
equation with a > 3/2, and showed that 

" exp[-u' 3 /(3Bi)] if u' oo 
w'~ 5 / 2 if u' —>■ — oo 

Clearly, the probability density of u' is nongaussian and has a power-law tail for negative u' . Boldyrev fl4|,|ll| extended 
Polyakov's method to the range l/2<er<3/2 and calculated P(u'). He derived that x = l + 2er/3for 1/2 < a < 3/4, 
and x — 3/2 for a > 3/4. Boldyrev |l4|] suggests that the above formula may also be valid for < a < 3/2. 

Gotoh and Kraichnan jlTJ] found the exponent of P{u') for the negative v! to be —3. Recently, E and Eijndcn 
jlTj and Kraichnan |TJ| argue the exponent to be —7/2. Gurarie and Migdal and Balkovsky et al. [^0| have 
applied instanton solutions for solving probability density for positive u' and obtained qualitative agreement with the 
numerical results reported by Yakhot and Cheklov pl| . In this paper we will compare our findings with the above 
mentioned results. 

The roughening exponents for KPZ equation xkpz and (3kpz have also been calculated by Zhang by replica method 
p2[ , and by Hentschel and Family using scaling arguments |^3| . They find a close agreement with the results of Medina 
et al. p0| . Numerically, Peng et al. |24] have calculated the exponents xkpz and (3k pz using finite-difference method, 
and Amar et al. |^5| and Meakin and Jullien |2(|^7j have calculated using lattice simulations. Their results are in 
qualitative agreement with the RG predictions. 

The energy cascade rate is one of the important quantities of interest in the statistical theory of turbulence. In this 
paper we have computed the energy flux of the Burgers equation for — 1 < a < 6. We find a constant flux for all a 
beyond 1/2. In our simulation we also varied the values of the parameters (v, A, D) and analyzed its effects. We find 
that there are interesting crossover from KPZ to Edward Wilkinson (EW) equation depending on the values of the 
parameter. These crossover results are reported in the Appendix A. 

The outline of the paper is as follows: We restate the structure function and P(u') of the noiseless Burgers and KPZ 
equations in section 2, and the energy spectrum of these equations in section 3. From the nature of the noise-noise 
correlation, we divide the a regime in three parts. This result is discussed in section 4. Section 5 and 6 contain the 
simulation method and results respectively. Here we calculate the energy spectrum and structure functions of the 
Burgers and KPZ equations for various values of a. The probability densities of u' and u is also discussed in this 
section. Section 7 contains a brief discussion on the cascade rate of the noisy Burgers equation. Section 8 contains 
discussion and conclusions. 



II. STRUCTURE FUNCTION CALCULATION FOR NOISELESS CASE 

We will show later that the Burgers equation with strongly correlated noise has behaviour similar to the noiseless 
Burgers equation. Therefore, we will discuss the noiseless equation briefly before embarking on the noisy Burgers 
equation. Burgers Q solved the Burgers equation exactly at large time t for vanishing viscosity. The velocity profile 
at a given time is linear except for the sharp discontinuities at the shock positions (see Fig. 1). The structure function 
of the noiseless Burgers equation S q (r) has been calculated using this solution (see Gotoh [7|, Bouchaud et al. Q and 
others). Assuming the ensemble average to be the same as the spatial average, Gotoh obtained 

1 f L 

Sq( r ) = ~f I \u(x + r) — u(x)\ q dx 
L Jo 

i 

where L is the length of the box, fii is the shock strength (defined as the velocity difference across the shock) of the 
ith shock, and £j is the position of the ith shock. The second term of the equation is due to the discontinuities at the 
shocks. For small r and q > 1, the second term will dominate the first one. Hence, for small r, S q (r) oc r^i 9 with 
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Cf™" 9 = 1. (2.2) 

Note that this relationship is valid for 5 <C r -C A, where 5 is the shock width, and A is the average distance between 
two shocks. For — ^ 0, (5 is finite but small (for details, see Saffman [p8| ). Since u is continuous with in the shock, 
S q (r) = r g for r -C 6. For large r, <S g (r) is not proportional to r because of (/ij — r/i) 9 term of Eq. (|2.l[) . 

A point is in order here. Tatsumi and Kida H showed that the number of shock fronts decrease with t as £~ 7 , 
where < 7 < 1. Hence in the asymptotic state, there will be only several shocks, and the distance between the 
shock fronts will be of the order of box size. Therefore, asymptotically A will of the order of the box size L, and r 
can be comparable to L. 

As mentioned in the introduction, structure function is closely related to intermittency. Since the exponent ^ ur 9 
is not proportional to q, the noiseless Burgers equation is classified as intermittent system. However, a point to 
note is that the velocity of the noiseless Burgers equation is not random. From the velocity profile it is clear that 
P(u) = const for u between u m in and u m ax- The slope u' is a constant (c) for all u except within the shock region 
where v! is large but negative. Therefore, P(u') will be a sum of a delta function at u' — c and a small spiky function 
at large but negative v! (see Fig. ||). Since u is not random, it is somewhat confusing to call the signal as intermittent. 
However, it is common practice to classify the noiseless Burgers equation as an intermittent system. 

For fluid turbulence it has been shown that the intermittent velocity field has a multifractal distribution. This was 
demonstrated by Meneveau and Srinivasan @,[30| using the cascade model. We will apply the same model to Burgers 
turbulence for a further understanding of the intermittent nature of the Burgers solution. 

In the cascade model of Meneveau and Sreenivasan |^9|,^0) for fluid turbulence, the nonlinear energy flux E r is 
distributed unequally between two smaller eddies. The flux E r at length scale r is divided into fractions p\E r and 
j)iE r to two smaller eddies of length r/2, and so on. The generalized-dimension of the resulting multifractal is given 
by 

D q =\og 2 {p\+piyl^. (2.3) 
The constants D q are related to the exponents of the structure function of the energy flux Ef. in eddy of size r by 

I>; = E1(I) <S - ,,D \ (2.4) 

where the sum is taken over all the eddies at the nth stage. When p\ ^ pi-, the flux function is unevenly distributed 
after several bifurcations of the original eddy, and that yields bursty behavior for E r . A closer inspection shows that 
the Burgers solution corresponds to p\ — > and P2 — > 1, which implies that D q rs 0. Hence, for Burgers turbulence, 
only one among the 2™ eddies has all the flux after n bifurcation of the original eddy. This corresponds to maximal 
intermittency or maximal multiscaling. Using Eq. (5.6) of Meneveau and Sreenivasan fjCfl , we get 

Cf^ = (g/3-l)D, /3 + l = l, (2.5) 

a result consistent with Eq. ( |2.2| ). Note that when pi = P2 = 1/2, E r is constant and D q = 1, and there is no 
intermittency. 

In Burgers equation, the dissipation occurs only at the shocks, which are narrow regions. So, as we traverse along a 
line, there are regions of no dissipation, then suddenly a short region (shock) with intense dissipation appears. Accord- 
ing to this observation, the noiseless Burgers equation exhibits strong intermittency, in fact maximal intermittency, 
in the light of Meneveau and Sreenivasan's model (2^j3^] . 

The structure function of the noiseless KPZ equation can be easily calculated using the relationship h(x, t) = 
— J x u(x' , t)dx' . The surface profile h(x) at a asymptotic time (shown in Fig. |l|) is 

Mx t) = { ~i + f(*-&) for 5(6-1 + &) < x < & 

\ -ft + 3! F{x-t i ) for & + 
After some algebra we can obtain the structure function T q (r) for the KPZ equation: 



1 , r y [(r? i+2 - 2&+i) 9+1 - (r?i+i - 2&+i) 94 j 
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The solution of the KPZ equation has cusps at positions where shocks appear in Burgers equation. These cusps 
yield contributions proportional to r q+1 to the structure function. Hence, to a leading order, T q (r) oc r q . Therefore, 
C!f PZ = Q- Note that smooth continuous curves also yield £ g = q because h(x + r) — h(x) « h'r for small r. Since the 
exponent Cg is proportional to q, the solution of the noiseless KPZ equation will be classified to be non- intermittent. 
However, note that the probability density of h! — —u is flat (nongaussion), hence indicating intermittent behavior 
for the noiseless KPZ equation. Clearly there appears to be a contradiction because the structure function indicates 
non-intermittency, but the P{h!) indicates intermittency. The resolution of this apparent contradiction is given below. 

The systems under investigation for intermittency have typically a random output. In non-intermittent systems, 
the probability density of the output is gaussian, i.e., probability of finding a large signal decreases as exp(— x 2 ). 
However, in the intermittent systems, the probability density of the random output signal deviates from the gaussian 
behaviour. The key is that the system under investigation for intermittency must have a random output, which is 
not the case for the noiseless KPZ and Burgers equation. In that sense, it is somewhat meaningless to ask whether 
the noiseless Burgers and KPZ equations exhibit intermittency or not. This is the reason why the conclusions from 
the structure function and those from the probability density differ for the noiseless KPZ equation. In literature, 
however, the noiseless Burgers equation is classified as an intermittent system. 

In the follo wing section we will derive the energy spectrum of the noiseless Burgers and KPZ equation using 
Eqs. (OB. 



III. ENERGY SPECTRUM OF NOISELESS BURGERS AND KPZ EQUATION 

We briefly discuss the energy spectrum of the noiseless Burgers and KPZ equation because they will be compared 
with the results obtained for the equations with coloured noise. We restate the earlier results by Kida Gotoh [[r|, 
and Bouchaud et al. |J. We can obtain the energy spectrum E u (k) =< \u(k)\ 2 /2 > using Ss(r). Clearly, 

< u(x + r)u(x) >=< u 2 > --S 2 (r). (3.1) 



Therefore, 



I f°° 

E(k) = — 1/2 < u(x + r)u(x) > exp(—ikr)dr 

L J-oo 



= S ^^- + ^^ + (Lk) (3.2) 

i 

Hence, the noiseless Burgers equation has energy spectrum E u (k) oc fc~ 2 for k > 0. In the above equation, the Fourier 
transform of r 2 is d 2 5k,o/dk 2 pjj |. 

A similar procedure using T^r) will yield E h {k) =< \h(k)\ 2 /2 >, 

^ = 4,0^+^4,0 + ^- (3.3) 

where B and C are constants |3l]]. Hence for k > 0, we obtain E h {k) oc fc -4 for the noiseless KPZ equation. 

Most papers in the past implicitly assume that if T^r) oc r 2x , then E h {k) should be proportional to fc~ 2x_1 . Clearly 
this does not hold for the noiseless KPZ equation (check: x = 1; but E h (k) oc fc~ 4 ). This apparent contradictions, 
which also occurs for a > 1, can be resolved using the following arguments. 

Suppose E(k) = Ak~ 2x ~ x . The second-order structure function S^ir) is given by 

I f L 

S 2 (r) = T dx < \u(x + r) - u{x)\ 2 > 
L Jo 

/>oo 

,2/ 



= 16 / dkE(k) sin 2 (kr/2) 
Jo 

= 16/ dkAk- 2x ~ 1 sm 2 {kr/2). (3.4) 
Jo 

We are interested in the leading order behaviour of S 2 {r) for small r. When X < 1, the above integral converges, and 
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S 2 (r) cx l6Ar 2x / dss' 2 ^ 1 sm 2 (s) 
Jo 

m ACir 2 ^, (3.5) 

where G\ is a dimensionless constant. When x > 1, the integral diverges from below (small k). However, we can cure 
this divergence by choosing the lower limit of the integral to be 2tt/L, which yields 

/>oo 

S 2 {r) cx 16A / dkk- 2x - 1 {kr/2) 2 

J2w/L 

« AC 2 (r/L) 2 L 2 *, (3.6) 

where C 2 is a constant. 

Hence, when x < 1, E{k) cx £; _2x_1 and S 2 (r) cx r 2x as expected, however when \ > 1, -E(fc) cx fc _2x_1 and 
S 2 (r) cx r 2 . This analytical results are seen in our numerical simulations to be described below. 

After the discussion on the noiseless equation, now we turn to the noisy Burgers and KPZ equations. 



IV. NOISY BURGERS EQUATION: VARIOUS RANGES OF PARAMETERS a 

In this section we will attempt to divide the parameter range of cr according to the properties of the solution. From 



Eq. (1.1) the noise spectrum |C(^)| i s proportional to k . Using this we can derive the noise-noise correlation as 



B S(r) for a = 

B 1 -Cir 2a - 1 for0<cr<l/2 
(((x,t){(x + r,t)) ~ { B 2 -C 2 \og{r) for a = 1/2 (4.1) 
B 3 -C 3 r 2a - 1 for 1/2 < a < 3/2 
B 4 - Cir 2 L 2a - 3 for a > 3/2 

where L is the length of the system, and Bi and Ci are constants. For — 1 < a < 0, |C(^)| 2 ~~ ¥ 00 f° r large k, hence, 
the inverse Fourier transform which yields < ((x, t)((x' , t) > is not defined in this regime. However, Kardar et al. |J 
and Medina et al. jl0| have solved the KPZ equation for this regime using \f(k)\ 2 which is well defined for large k. 
We call the region — 1 < a < as (a). 

It can be deduced from the above discussion that the noise-noise correlation increases with the increase of a till 
fj = 3/2. Beyond a = 3/2, the correlation is proportional to 1 — Cr 2 for all a. Therefore, it is expected that the 
h — h and u — u correlation would increase with the increase of a till 3/2, beyond which the behaviour is expected to 
be somewhat similar. Keeping this in mind we have divided the a range beyond in two regions: (b) < a < 3/2, 
and (c) a > 3/2. Well defined shocks develop in the parameter range (c) due to the large noise-noise correlation, and 
these shocks determine the exponents. In the intermediate parameter range (b), the exponents change slowly from 
RG dominated values to shock dominated values. The details are given in the next section. 

Polyako v |2|| h as analytically solved Burgers equation with cr = 3/2 using the methods of quantum field theory. 
Boldyrev jl4j|l5| extended Polyakov's method to the range 1/2 < a < 3/2 and derived the spectral exponents and 
the probability densities. Boldyrev derived that x — 1 + 2<r/3 for 1/2 < cr < 3/4, and \ = 3/2 for a > 3/4. Boldyrev 
jlij suggests that the above formula may also be valid for < cr < 1/2. In the following section we will compare our 
numerical results with the theoretical predictions of Medina et al. Polyakov and Boldyrev |l4| , |l5| ]. 



V. SIMULATION METHOD 



Our calculations in this paper have been done using direct numerical simulations based on pseudo-spectral method. 
This method, commonly used in turbulence simulations, is expected to perform better than finite difference scheme 
because the derivatives can be calculated exactly in the spectral method |$2| . The finite difference scheme was adopted 
by Moser et al. |33| and Peng et al. |24j in their simulation of KPZ equation for — 1 < a < 0. 

We solve the KPZ equation in one dimension. The details of the simulation are as follows. A box of size 27r is 
discretized into N — 1024 divisions. The KPZ equation is solved in Fourier space. However, to compute the nonlinear 
term, we go to real space, perform multiplication, then again come back to Fourier space. We time advance the 
Fourier components h(k,t) using Adam Bashforth time marching procedure with flat surface as an initial condition. 
Two-third rule is used to remove aliasing |32|. For details of the simulation refer to Canuto et al. |32j and Verma et al. 
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In our simulation we also add hyperviscosity term («V h) to the KPZ equation to damp the large wavenumber 



modes strongly. The hyperviscosity term does not affect the intermediate scales which is of our interest 34 1; this is 
because the V 4 /i and higher order derivative terms are irrelevant in the renormalization group sense (cf. Barabasi 
and Stanley jl2| and references therein). 

In our simulation we take v and k to be very small. Note that large v correspond to Edward- Wilkinson (EW) 
equation. The dimcnsionless parameters used in our simulations are 

A = 1.0, 

v = W- 5 , 
k = 10~ 6 , 
(2n) 2 D = 10~ 3 , 

^= 1/2000 (5.1) 

with one exception. For a = — 1, we choose A = 0.1 for the stability of the code. In the Appendix A we have varied 
values of the parameters v and D and shown crossover from KPZ behaviour to EW equation. 

The values of a used in our simulation are -1, -0.85, -0.75, -0.60, -0.50, -0.25, 0, 1.25, 1.50, 1.75, 1, 2, and 6. 
We have time evolved the equation till 15 nondimensional time units. We find that the system reaches saturation 
in approximately 8 to 10 time units. For reference, It: time unit corresponds to one eddy turnover time in fluid 
turbulence. For ensemble averaging, we have performed averages over 100 samples which start with different random 
seeds for the noise. We have used rani of numerical recipes p5[ as our random number generator. Each computer 
run for 15 time units and 100 samples takes approximately 7 hours on a Pentium machine (150 MHz). 

We have calculated \ an d P using the simulation data. The ensemble average < . > have been obtained by taking 
averages over 100 runs. The width W(L, t) = [< h 2 {x)/N >] 1 / 2 grows as a power law in time, i.e., W(L.t) cx t^ KPZ , 
in the early stages of growth. We obtain (3kpz by fitting a straight line in log — log plot of W(L, t) vs. t over a range 
of t = 0.2 : 2.5. The other exponent xkpz is obtained from the asymptotic < \h(k,t)\ 2 > averaged over 100 runs. 
We find that powerlaw holds (|/i(fc)| 2 cx k~ 2xKPZ ^ 1 ) for k = 10 : 8; this range corresponds to the inertial range of 
turbulence. We perform the averaging at steady-state (t — 15). In Fig. || we plot < |/i(fc)| 2 > vs. k for a = —1, 0, 1/2, 
and 2. The lines of best fit are also shown in Figure. The computed values of the exponents are listed in Table 1. 
Our estimate of the error in the exponent is roughly 0.05. 

We have also calculated the structure functions S q (r) and T q (r) using the simulation data. The variable r ranges 
from 1 to N/2. Near the boundaries we calculate Am or Ah using wrap around scheme, i.e., x + r is taken as 
mod{x + r,N). we have reported the structure function exponents at the steady-state (t = 15). 

The probability densities P(h), P(u), and P{u') have been calculated by averaging the histograms over 100 runs 
at every 0.2 time interval from initial time of 5 units to the final time of 15 units. Even though the solution has not 
reached the steady state at t = 5, for more sampling we have taken the time interval from t — 5 to t = 15. 

In the following section we will describe the results of our simulation for various degrees of noise correlations, i.e., 
for different as. 



VI. SIMULATION RESULTS 



As discussed in section IV, we divide the range of parameter a in three regions: (a) — 1 < a < 0, (b) < a < 3/2, 
and (c) a > 3/2. Fig. |i| shows numerical u(x,t) for various us. We find that noise dominates when a is negative. 
However, shock-like structures are clearly visible for a — 2,6 (compare with Fig. [I]). In the the intermediate range 
of a (region (b)), both noise and shock structure coexist. We show below that the structures present in the profile 
contribute significantly to the determination of the spectral and intermittency exponents. 

Regarding energy spectrum, our numerical results are in good agreement with the RG predictions for — 1 < er < 
or < p < 1. There is a gradual transition from xkpz = 1 to xkpz = 3/2 as a increases from to 3/2. For a > 3/2 
we find that xkpz = 3/2. 

Regarding probability density, we find that the probability densities P(h) and P(u) are gaussian for all er's. Fig. || 
shows the plot of P(u) vs. u along with the best fit gaussian curves. Clearly P(u) is gaussian for all the ers. Therefore, 
the exponents of the structure functions T q (r) are expected to be proportional to q. 

Fig. ^| shows a typical log-log plots of S q (r) and T q (r) vs. r (here shown for a = 6). The shock region is r < 10 
where both S q (r) and T q (r) are proportional to r q because u(x) and h(x) are continuous here. The region of our 
interest is r = 10 : 80 where both S q (r) and T q (r) are power laws. The exponents for this range of r, ^ ur 9 anc | (^ pz , 
are listed in Tables 1 and 2 respectively. Fig. [7| shows C!f PZ vs - <Z plot, and Fig. || shows ^ ur s V s. q plots. Clearly, 
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C^ PZ is approximately proportional to q as predicted in the previous paragraph. Hence, our results regarding P(u) 
and T q (r) are consistent. 

We find, however, that P{v!) deviates significantly from the gaussian behaviour as we increase a from 0, thus 
signalling an intermittent behaviour for u. The details of our results for various ranges are given below. 



A. -1 < a < 



In Table 1 we have listed our numerical values of xkpz and (3kpz- For comparison we also list the predicted values 
of Medina et al. |nj below 



1 for -1< a < - ' 

. for — — 

3 4 



1 + 1(7 for -f < cr < 4 



and 



for -1< a < 



3 



P=\ h*> Zl^n 4 (6-2) 



3 _ 2CT for-f<(7<0 



Medina et al.'s results are based on RG scheme that breaks down beyond a = 0. 

We compute xkpz by a straight line fit to the log — log plot of < \h(k)\ 2 > over k — 10 : 80 at t = 15 (see Fig. ||) 
Note that \h(k)\ 2 oc k~ 2x ~ x . The xkpz$ listed in Table 1 show that that our numerical xkpzs are in close agreement 



with Medina et al.'s [|10| and Chattopadhyay and Bhattacharjee's 11 theoretical predictions. Our findings are also 
consistent with earlier simulation results by Amar et al. p5| , Peng et al. and Meakin and Jullien |2^,^]. For 
cr = —1/4, Hayot and Jayaprakash |3| have reported a crossover in the wave number space — from RG dominated 
region for small k to free field behaviour for large k. We do not find any such crossover in our simulation. Regarding 
the Burgers equations, XBurg = Xkpz - 1- 

Regarding (3kpz calculations, our results are agreement with Medina et al.'s predictions for a = — 1 : —3/4. 
However, our exponents differ significantly with Medina et al.'s exponents for higher as. For example, for a = 0, we 
obtain (3k pz = 0.45 contrary to the predicted (3k pz = 1- The reason for this discrepancy is not clear to us at this 
point. Regarding Psurg, we find that it is approximately 0.5, consistent with Hayot and Jayaprakash's findings 

The probability density P(u') vs. u' for a = —1, —1/2, is shown in Fig. [)] on a semilog plot. The figure shows 
that P(u') is gaussian for a — — 1, but it deviates from the gaussian behaviour with the increase of a; the deviations 
becomes more and more prominent for higher a. At a — 0, P(u') ~ 0.01 * exp(— 0.14|m'|) for negative u' (the curve 
of the best fit shown in the Fig. [|). Beyond a = 0, a powerlaw tails start appearing for negative u' . 

Table 2 and Fig. || show that ^ urg is close to zero for a < 0. The energy spectrum < \u(k)\ 2 /2 >oc k~ 2xKPZ+1 . 
Therefore, using the arguments of section 3 we can easily show that < \u{x + r) — u{x)\ 2 >cx r 2i ° iKPZ ~ 1 ^ for a < 0. 
For a — 0, it can be easily shown that < \u(x + r) — u(x)\ 2 >cx log(r). We find in our simulation that (2 PZ i s 
positive contrary to the above prediction. This deviation may be because of intermittency (due to the presence of 
small shocks). It is also possible that at large r (r i=s L), the other sub-critical terms may become comparable to the 
leading order term and may change the exponent. Quantitative calculation of the exponents for this range of a is 
beyond the scope of this paper. 



B. a > 3/2 



One of the important aspect of this paper is the discussion of the energy spectrum, probability density and structure 
functions of the Burgers equation with large a. As shown in the Fig. ^[ shocks are prominent for this range of a. 
Due to these shocks, the exponent xkpz ~ 3/2 for cr > 3/2 (see Table 1 and Fig. ||). Note that xkpz = 3/2 and 
XBurg = 1/2 for the noiseless Burgers equation due to the presence of shocks. 

The existence of shocks for large a can be argued from the noise-noise correlation discussed in section IV. For large 
cr, there is long-ranged noise-noise correlation. Physically, large fluid parcels are moved around by this noise. As a 
consequence, there will be regions where the parcels forced by appositely directed noise collide with each other and 
create "strong" shocks. Hence, it is not surprising that strong shocks will be generated for large a. These shocks 
determine the dynamics and the spectral indices of Burgers and KPZ equations. Therefore, the noisy KPZ equation 
with large a yields the same energy spectrum and multiscaling exponents as the noiseless KPZ equation. 

The noise- noise correlation is proportional to 1 — Cr 2 for all a > 3/2, hence the spectral and multiscaling exponents 
are expected to be somewhat similar for all the a beyond 3/2. This is borne out by our numerical simulation. This 
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is consistent with the trivial observation that for large a only k = 1 mode is effective. It is interesting to note that 
k = 1 mode (i.e., C( x ) = sin(x)) also yields noise- noise correlation proportional to 1 — Cr 2 , and approximately the 
same exponents as those with a > 3/2. 

We have also analyzed the stability of the shocks at a preliminary level. For large a, the significant contribution to 
the dynamics comes only from the k — 1 mode of the noise. Since the first mode (k = 1) is noisy, the the formation 
of a single shock delayed for some time because of the movement of the zero of the noise signal. However, after the 
shock is formed, it shifts around by only a small amount because the impulse due to the random noise is not strong 
enough to move the shock by a large distance. 

We have also calculated (3 for both Burgers and KPZ equations when a > 3/2 and found them to be approximately 
1/2 . This result is consistent with findings of Hayot and Jayaprakash 

Regarding the intermittency exponents, as shown in Fig. |g, Q^ ur 9 ~ \ fo r all q for a > 3/2. This result is consistent 
with existence of shocks (see Fig. [| and section II). For large r, S q (r) is not proportional to r. This is because of the 
term similar to (/i^ — r/t) q of Eq. ( |2.1| ), which becomes important at large r. The effects of this subcritical term is 
evident in the intermittency exponents listed in Table 2. We had anticipated ^ ur 3 to be 1, but we consistently find 
them to be less that 1. The same trends were observed for the noiseless case. 

Fig. shows the probability density P{u') for various crs. The curves of best fit are also shown in the Figure. 
The probability density is clearly nongaussian for all crs shown in the Figure. In our simulations we find that for 
all (7 > 3/2, P(u') oc u'~ a (a — 2.2 — 2.5) fits quite well with the numerical P(u') for the intermediate range of the 
negative u'. This result is in agreement with the Polyakov's theoretical prediction for a > 3/2 J2|. Curiously, P(u') 
for v! > is also nongaussian, but it does not match with the Polyakov's predictions for positive v! . 

We have shown in this section that both noiseless Burgers equation and the Burgers equation with strongly correlated 
noise have SJr) cx r. However, the probability densities P(u) and P(u') are very different for these two cases (see 
Figs. H and [h]). This difference is because of the different nature of the noise. 

The entries in Table 2 show that C^ pz = q. This can be argued as follows. Since xkpz ~ 3/2, the second order 
structure function T2(r) will be proportional to r 2 , or C,2 PZ — 2 (see Section III). This is seen in our simulation for 
a > 3/2. The Fig. [| also shows that P(u) is gaussian. Therefore, T q (r) cx \T2(r)] q l 2 cx r q , hence C^ pz = <?> a result 
seen in our numerics. Thus, structure function and energy spectrum calculations are consistent. 

After discussing the shock dominated region of a > 3/2, we turn to the region < o < 3/2. 



C. < a < 3/2 

The Xkpz calculated by our numerical simulation for the range < cr < 3/2 is listed in Table 1. The exponent 
increases from 1 and saturate at 3/2. That is, there is a gradual shift from RG dominated exponents to shock 
dominated exponents as we vary cr. The profile u(x,t) of Fig. 4 shows that both fluctuations and shocks coexist in 
this range, with shocks becoming more and more important as a increases. 

Boldyrev 14 |l^] has calculated xkpz using methods of quantum field theory and derived that xkpz = 1 + 2cr/3 



for 1/2 < a < 3/4, and xkpz = 3/2 for a > 3/4. Boldyrev |14| suggests that the above formula may be valid for 
< cr < 3/2. Comparison of Boldyrev's predictions with the entries of Table 1 shows that the predictions work quite 
well till < a < 1/2. After that there is a significant deviation, and the exponents appear to be outside numerical 
error bars. We need to probe the region 1/2 < cr < 3/4 carefully to reach a definite conclusion. 

We presume that both fluctuations and the embedded structures are important in the determination of the spectral 
index. One would need to combine contributions from the structures and fluctuations to obtain the energy spectrum 
p6f . In Appendix B we sketch an elementary framework when both structure and fluctuations are present in the 
system. 

The Burgers equation for cr = 1/2 has been studied extensively by Chekhlov and Yakhot Q. They find in their 
high resolution simulation that |u(fc)| 2 cx fc -5 / 3 . The spectrum in our low resolution simulation is |u(fc)| 2 cx fc -154 
{xkpz — 1-27), which is close to Chekhlov and Yakhot's [Q result. Chekhlov and Yakhot have argued for a 
constant cascade of energy in the wavenumber space in this case and claimed that the Kolmogorov-like energy spectra 
is due to the constancy of the energy flux. However, we find the cascade rate to be constant for all a > 1/2, but the 
spectral index of the energy spectrum is in the range of 5/3 to 2 (to be discussed in Section 7) . Hence, the argument 
that the constant flux yields Kolmogorov's spectrum is incorrect. 

Regarding the probability distribution, in Fig. |l^ we plot P(u') vs. v! along with the curves of best fit. It is found 
that P(u') is a powerlaw for v! < 0. For a = 1.25, P(u') cx u'~ 2 - 86 , but for a = 1/2, P(u') w w'~ 2 - 2 , with an error 
of approximately 10% for both the cases. Curiously, Polyakov's predictions (cr = 3/2) for positive u' appears to be 
applicable only for cr = 1/2; here P(u') cx exp(— Cx 3 ). 
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The roughening exponents xkpz lies between 1 and 3/2. Therefore, C^ FZ should be proportional to q as argued 
in the previous subsection. Our numerical simulations yield approximately the same exponents (see Table 3). Hence 
our energy spectrum calculations and the structure function calculations are consistent. 



D. Summary of simulation results 



From the above discussion we see that Burgers and KPZ equations are well understood for —1 < a < and a > 3/2. 
In the intermediate range < a < 3/2, Boldyrev's predictions appear to explain the numerical data in part of the 
regime; further analytic and numerical work is required in this regime. 

For all a > —1/2, xkpz + z (z = xkpz / Pkpz) deviates significantly from 2. In fact, for larger er, Xkpz + z — 
3/2 + (3/2)/(l/2) = 9/2, quite far from 2. Note that Medina et al. |1Q[ argue that the identity xkpz + z = 2 is 
a consequence of Galilean invariance. However, Hayot and Jayprakashjsj] , Polyakov Q and others have speculated 
violation of this identity due to the presence and motion of the shocks. Meakin and Jullien's |26 2?|] results also violate 
Xkpz + z — 2 condition for a range of a. They find that for a = —1/4, xkpz + z = 2.37, quite different than 2. The 
violation of Galilean invariance even in the region — 1 < a < signals importance of structures 

We also show that P(u) is gaussian for all a. However, P(u') is gaussian at a = — 1, but continuously changes to 
exp(— x) then to a powerlaw behaviour. Therefore, u field of Burgers equation exhibits intermittency. The degree of 
intermittency depends on a. 

The energy flux play an in important role in turbulence analysis. In the following section we briefly report energy 
flux studies for the Burgers equation. 



VII. ENERGY FLUX OF NOISY BURGERS EQUATION 

We derive an equation which gives us the energy transfer from the region |fc| < K. The derivation of the energy 
equation from Eq. (1.1) and averaging yields |37J 



— J E(k)dk = - J 2vk 2 E(k)dk 



K 



ft < u*{k) 



f K 

>+ SR< u*(k)f(k) > 
k Jo 



From this equation, clearly the energy dissipation in the wavenumber sphere of radius K is 

r K 

D K = 2vk 2 E(k)dk, 
Jo 

and the energy flux coming out the wavenumber sphere of radius K is 

-JC r / g 



Hi 



ft < u*(k) 



o 



ET I -.72 



> . 



(7.1) 



(7.2) 



(7.3) 



It can easily seen that for the random noise, the energy supplied by the forcing to the wavenumber sphere of radius 
KisM 



F K 



K 



ft < u*{k)f{k) > 



K 



< \f(k)\ 2 > ■ 



From the Eq. ( |7.1| ) it is clear that at the steady state 

F K = U K + D K . 



(7.4) 



(7.5) 



We have plotted Dk,Hk, arid Fk for various a. Fig. [Tl] shows the plots for a — 0. Here, we find that Dk ~ Fk 
till k ph 250, but Hk <C Fk except for small k. Similar results are obtained for all — 1 < a < 0. 

Fig. 12 shows the plots of Dk,Hk, and Fk for a = 1/2. The forcing rate Fk is proportional to \og(K). We also 
find a range of K for which II^ is constant. These results are consistent with the findings of Chekhlov and Yakhot [Q. 
Fig. 03 shows the plot of the above quantities for a — 6. Here again II^ is constant for a range of K. An interesting 
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point to note is that for a — 6, Fx is constant beyond K — 10 or so; this is because effectively only first few modes 
are forced when a is large. Similar results are obtained for all a > 2. 

We find that the energy flux is constant for the noisy Burgers equation for all a > 1/2. However, the energy 
spectrum varies from fc -5 / 3 to fc~ 2 as we increase a from 1/2 to 3/2 and beyond. Hence our results show that 
constancy of energy flux is not a sufficient condition for the Kolmogorov's energy spectrum in noisy Burger equation. 
This indicates that Chckhlov and Yakhot's jj| argument that the constant energy cascade rate for a = 1/2 implies 
Kolmogorov's energy spectrum is incomplete. 

We find that for a > 2, the flux rate Hk ~ 2 x 1CP 3 . This numerical value is consistent with value obtained using 
the formula derived by Saffman p8[ 

IIk « 4 (7.6) 

where \x is the velocity jump across the shock. In our simulations, fi s» 0.3 for large a. Also, the Eq. (|7.4j) yields the 
same numerical value {2D — 2 x 10~ 3 ). 



VIII. DISCUSSION AND CONCLUSIONS 



In this paper we have numerically calculated the energy spectrum, structure function, and probability densities 
(P(h), P(u), P(u')) for the KPZ and Burgers equations in the presence of correlated noise. The Burgers equation 
with a strong correlated noise (a > 3/2) has distinct shock structures. Given this we have trivially worked out the 
roughening exponents for this range: XBurg = 1/2, XKPZ = 3/2. The \ calculated from the numerical simulation 
matches quite well with this result. It has been theoretically shown by many researchers that Cf vrg « l when shocks 
are present. Numerically we find the same Cq Vrg f° r a > 3/2. The probability density P(u') oc u'~ a with a — 2.5±0.3 
for negative u'; this result is in agreement with the Polyakov's theoretical predictions. However, for positive u', our 
results do not match with Polyakov's predictions for this range of a. Note that our results differ from those of E and 
Eijnden jL7| and Kraichnan jig] who obtained a = 7/2, and Gotoh and Kraichnan ]^6) who argue for a = 3 

In the a regime —1 < a < 0, our numerical xkpz matches with the RG predictions. The exponent Pkpz, however, 
saturates at 0.5 contrary to the RG predictions. For < a < 3/2, the exponent xkpz varies smoothly from 1 to 3/2. 
Boldyrev's predictions appear to match with our numerical results for < a < 1/2. After that there is a significant 
deviation. Our claims, however, are not on a very strong ground because of large uncertainties, and hence, more 
exhaustive simulations are needed to reach definite conclusions. 

We find in our simulation that P(u) is gaussian for all as and also (J^ FZ oc q. Therefore, according to the convention 
discussed in the introduction, the signal h will be termed as non-intermittent for all as. Note, however, that according 
to the same convention, u is as intermittent for a > 1/2. These two statements appear to be contradictory because 
u and h are related by u = — V/i. The apparent contradiction is quite simple to resolve. The velocity signal u, being 
the derivative of h, is much more singular than h because of the presence of the shocks. That is why probability 
density P(u') is nongaussian (intermittent u), even though P(u) is nongaussian (non-intermittent h). This observation 
cautions us to choose a right variable while investigating the system for intermittency. Note that for fluid turbulence 
also, P(u) is gaussian, but P(u') is not, and the variables used for structure function is u (there is no corresponding 
h anyway). 

The nonexistence of intermittency in KPZ equation will have relevance to other models of surface growth. Sneppen 



and Jensen [J38| and Tang and Leschhorn | 
of quenched disorder. They find no spatia 



39j calculated the structure function in surface growth equation in presence 
multiscaling in their system; this result could be related to the conclusions 
described in our paper. Sneppen and Jensen pS| ], and Tang and Leschhorn p^ ] find temporal intermittency in their 
model. We believe temporal structure function for the noisy KPZ and Burgers equation are important and will shed 
further insights into the dynamics of these systems. 

Krug jl0| and Kundagrami et al.lE) have investigated the existence of intermittency in the surface growth model of 
Das Sarma and Tamborenea (DT) ||42|j model and found multiscaling in it. Krug, however, finds absence of multiscaling 
in a variation of the DT model that was tilt independent. The DT model is related to linear Lengevin equation which 
is not expected to show intermittency, but Krug has attempted to relate intermittency in the DT model with the 
existence of relevant variables of renormalization groups. The connection of Krug's result to KPZ equation is not 
clear to us at this stage. Also, it is not clear whether the multiscaling in the DT model is an artifact of lattice effects, 
or it is due to appearance of relevant of nonlinearities in the continuum equation p7j| . Another important point to 
note is that u of fluid turbulence should probably be mapped to h', not h as done in Krug's and Kundagrami et 
al.'s |4lJ] papers. 
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In the regime < a < 3/2, both structures and noise coexist. Recently, Polyakov g], Boldyrev [ [14p5| ], Gurarie and 
Migdal 19 1, E and Eijnden |17|] , and Kraichnan (ll| have attempted to analytically solve for the spectral exponents 
and the probability density P(u') for this regime. Many of the recent attempts are based on methods of quantum field 
theory. Considering the important role played by the shock structures in determination of the roughening exponents, 
we believe that a calculation which incorporates both structures and noise will be very useful for such problems. 
The roughening exponent will get contributions both from fluctuations and embedded structures (see Appendix B). 
However, a careful analysis is required to isolate the individual contributions. Krishnamurthy and Barma ]36|] have 
isolated a moving pattern in the surface growth phenomena in the presence of quenched disorder. Their method may 
be applied here to separate the fluctuations from the structures. 

Bouchaud et al. [|S| have calculated the structure function of the Burgers equation in higher dimensions using the 
connection of KPZ equation to directed polymers. They showed that S q (r) ot r. In Appendix C we argue that using 
the shock structures of the higher dimensional Burgers equation, one can obtain S q (r) cx r. 

In this paper we have demonstrated the usefulness of structures in calculating the dynamical exponents of the 
system. The role of structures in dynamics is being studied in fluid turbulence, intermittency, self organized criticality 
etc. For example, in fluid turbulence Hatakeyema and Kambe {S| have used the vortex structures to calculate the 
scaling exponents. Therefore, discovery of the connections between the structures, fluctuations, and dynamics will 
yield interesting insights in the nonequilibrium phenomena around us. 
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APPENDIX A: EFFECTS OF PARAMETERS v AND D 



The roughening exponents of KPZ equation are usually stated without reference to the range of the parameter v 
and D. Usually it is assumed that v is small. This is in the same spirit as in fluid turbulence. While performing our 
simulations with correlated noise, we, however, found interesting changeover in the behaviour of KPZ even when v 
was small. We describe our findings below. 



As discussed in the main text, for the parameter described in Eq. (5T), we get the exponents shown in Table 1. 
However, in one of the test runs we fixed v at a somewhat higher value v = 0.05, and chose D = 0.001. For these 
parameters we found a transition for the energy spectrum from k~ 2p ~ 2 to k~ 2p and finally to k~ XKPZ , where xkpz is 
given in Table 1. To understand this transition, we need to look at the energy spectrum of Edward- Wilkinson (EW) 
and Random deposition (RD) models. 

EW equation is given by 

^M^V 2 Mx,i) + /0M). (Al) 
The above equation can be solved easily in Fourier space, which yields 

= (A2) 



Using 

(/(*, w')) = 2D k ^ (A3) 

l(L) — LO j 

we can easily show that 

(h(k, t)h* (k, t)) = k- 2p - 2 (A + B cxp(-2^fc 2 i) (A4) 
Therefore, EW equation at large t yields |/i(fc)| 2 oc k~ 2p ~ 2 . 
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If we substitute v = in the Eq. (Al), we obtain the Random deposition (RD) model [Q which is 



^ =/(*,*). (A5) 

Using a similar procedure as above, we obtain |ft.(fc)| 2 oc k~ 2p for this equation. 

Now we can explain the crossover from EW to RD and then to KPZ behaviour. In the initial phase, the width 
< h 2 > is small, hence the viscous term dominates both noise and nonlinear term, consequently EW behaviour is 
seen. As time progress and the width increases, the noise term becomes dominant and RD behaviour is observed. 
In this regime the width increases linearly in time fl^| . At a later time when the width become large enough, the 
nonlinear term takes over, and KPZ like behiour is observed. This is the reason why the spectral index varies from 
EW to RD, and finally to KPZ regime. 

Our findings regarding the crossover from EW to KPZ etc. show that we should be careful in the choice of 
parameters while simulating KPZ and Burgers equations. 



APPENDIX B: EXPONENTS WHEN BOTH STRUCTURES AND NOISE ARE PRESENT 



A surface profile h(x) can be split into two parts: 

h(x) = h s (x) + h f (x), (Bl) 

where hf(x) denotes the fluctuation, and h s (x) denotes the the embedded structure. By definition, < hf(x) >= 
where < . > denotes the ensemble average. Therefore, < h(x) >— h a (x). 

We can show that the width squared W 2 , which is obtained by taking the spatial average, is 



W 2 = ((h(x)-h) 2 ) (B2) 

= w 2 + w 2 

Hence, W 2 is the sum of the contributions from the fluctuating part and the structure part. Similarly, the second-order 
structure function T^r) is the sum of two contributions, 



T{r) = (\h{x + r)-h{x)\ 2 ) (B3) 
= T s (r)+T f (r). 

Hence, in a profile with a embedded structure, the exponent is determined by both structure and fluctuations. 

There is no structure in KPZ profile for a = —1, hence, the roughening exponent is totally determined the fluc- 
tuations. However, the profile for the noiseless Burgers equation has only shock structure, which determines the 
roughening exponent. In KPZ equation with correlated noise, especially for < a < 3/2, we have both structure 
and fluctuations. Therefore, we will have to do a careful analysis combining both of them to obtain the roughening 
exponents for these cases. We expect that similar analysis have to be carried out for the surface growth profile with 
quenched disorder that has both structure and fluctuations. 



APPENDIX C: STRUCTURE FUNCTION FOR BURGERS EQUATION IN D > 2 

The arguments of Gotoh Q for the computation of structure function of the noiseless Burgers equation can be 
generalized to higher dimensions. In two dimensions the projections of the shocks will appear as shown in Fig. |l4| . 
Towards the ends (labelled as Ai and Bi in Fig. |l4|) the shock strength tapers off and goes to zero. The velocity 
in between the shocks is a smooth function, but the velocity is discontinuous at the shocks. Therefore, as argued in 
section II, 

S q (r) = J |u(x + r) -u(x)| 9 dx 

KCi \G) q A+ c ^T,(^ q ( L ^ ( C1 ) 
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where A is the area of the system, \ii and Li are the shock strength and length of the «th shock respectively, and 
C\ , C2 are constants. The first term is obtained from the integral over all the space except the boxed regions near the 
shocks, while the second term is due to the integral in the boxed region. Essentially, we obtain |u(x + r) — u(x)| cx/i,; 
in the boxed region when x + r is in one side of the shock, while x is on the other side of the shock. This result yields 
S q (r) oc r for q > 1. These arguments for two dimensions can be easily generalized to three and higher dimensions 
(in three dimensions, the shock regions will appear as disks). Hence, in any arbitrary dimensions, noiseless Burgers 
equation has 

S q (r) cx r. (C2) 

Bouchaud et al. || have obtained this result in D dimensions by relating to KPZ equation to directed polymer; 
their result becomes exact as D — > 00. The strong intcrmittency of the Burgers equation is due to the large scale 
singularities of the velocity field at the shock regions, which are concentrated in D — 1 dimension (codimension 1). It 
was pointed out by one of the referee that the shocks could be fractal objects, whose fractal dimension Df could be 
greater than D — 1. The calculation of structure function for these objects is beyond the scope of this paper. 
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Figure Captions 



Fig. 1 The solutions of noiseless Burgers and KPZ equations. 

Fig. 2 The probability density P(u') vs. u' /u' max for the noiseless Burgers equation. 

Fig. 3 The energy spectrum (1/2) < |ft.(fc)| 2 > vs. k for the KPZ equation when a = —1,0, 1/2,2. The lines of best 
fit for the intermediate range of k are also shown. 

Fig. 4 The normalized velocity profile u(x) vs. x for Burgers equation when a = —1, —1/2, 0, 1/2, 1, 2, 6. The profile 
is noisy for — 1 < a < 0, has only well developed shocks for a > 3/2, and has both fluctuations and shocks for the 
intermediate range of a. 

Fig. 5 The probability density P(u) vs. u/u max for a — —1,0, 1/2 and 6. The best fit gaussian curves are also 
shown in the figure. 

Fig. 6 The structure functions S q (r),T q (r) vs. r for a = 6. 
Fig. 7 The exponent C q PZ vs. q for a = -1, -1/2, 0, 1/2, 2, 6. 
Fig. 8 The exponent C q urg vs. q for a = -1,-1/2,0,1/2,2,6. 

Fig. 9 The probability density P(u') vs. u'/u' max for a = -1,-1/2,0 and 1/4. Function 0.01 cxp(-0.14>'|) fits 
well for a = 0. 

Fig. 10 The probability density P(v!) vs. u'/u' max for a = 1/4,1/2,6. The best fit curves for negative u' are: 
M /-2.85 f or a _ u '~ 2 - 2 for a — 1/2, and u'~ 2 - 5 for a = 6. For positive v! , exp(— Cx 3 ) fits well for a = 1/2. 

Fig. 11 The Dissipation rate Dk , the energy flux Hk, and the forcing rates Fk vs. K for a = 0. 

Fig. 12 The Dissipation rate Dk, the energy flux , and the forcing rates Fk vs. K for a = 1/2. 

Fig. 13 The Dissipation rate Dk, the energy flux n^, and the forcing rates Fk vs. K for a = 6. 

Fig. 14 The projection of shocks of two-dimensional Burgers equation. In the middle, the velocity difference across 
the shock is large, while at the ends Ai and Bi, the difference is small. 
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TABLE I. The exponents a and /3 for KPZ equation for various cr's. The Medina et al.'s exponents are also listed for 
comparison. 
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TABLE II. The structure function exponents Cf urg 


of < \u(x + r) 


- u(x)\ q > for the Burgers equation for various cr's 
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TABLE III. The structure function exponents Cf PZ of < \h(x + r 


)-h(x)\ 


9 > for the KPZ equation 
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